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Abstract 

Abstract: High temperature expansions for the free energy, the susceptibil- 
ity and the second correlation moment of the classical TV-vector model [also 
known as the 0(N) symmetric classical spin Heisenberg model or as the lat- 
tice O(N) nonlinear sigma model] on the simple cubic and the body centered 
cubic lattices are extended to order /5 21 for arbitrary N. The series for the 
second field derivative of the susceptibility is extended to order /3 17 . We re- 
port here on the analysis of the newly computed series for the susceptibility 
and the (second moment) correlation length which yields updated estimates 
of the critical parameters for various values of the spin dimensionality N, 
including iV = [the self-avoiding walk model], N = 1 [the Ising spin 1/2 
model], N = 2 [the XY model], N = 3 [the classical Heisenberg model]. For 
all values of N we confirm a good agreement with the present renormalization 
group estimates. A study of the series for the other observables will appear 
in a forthcoming paper. 

PACS numbers: 05.50+q, ll.15.Ha, 64.60.Cn, 75.10.Hk 
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I. INTRODUCTION 



The continuing interest in the high temperature (HT) expansions for the statistical me- 
chanics of lattice spin models (or equivalently in the strong-coupling expansions for Euclidean 
lattice field theories) has in the last few years been a strong incentive for a substantial ex- 
tension of series in a variety of models. This valuable computational advance has been made 
possible not only by the large improvements of the computers performance and the rapid 
growth of their memory capacity in the last decade, but mainly by a more careful reconsid- 
eration of well known expansion techniques, which have yet to be fully exploited, and by a 
greater effort in devising and implementing faster algorithms . 

We have devoted this paper to the widely studied N— vector model on three- 
dimensional lattices, which is also known as the Heisenberg classical O(N) spin model or, 
in field theoretic language, as the lattice O(N) nonlinear sigma model. We recall that a 
HT expansion of the susceptibility \ to order (3 23 has been recently obtained || on the 
simple cubic (sc) lattice for N = [the self-avoiding walk (SAW) model 0] by a direct walk 
counting technique which cannot be extended to different values of N. In the N = 1 case 
[the Ising spin 1/2 model] the series 0{f3 21 ) on the body centered cubic (bcc) lattice for both 
X and the second moment of the correlation function /i2, obtained |||7j in the pioneering 
work by B. G. Nickel, remain still unsurpassed, but for iV > 1 the series now available are 
significantly shorter. The general situation before our work of which a brief, partial and 
preliminary account has already appeared in Ref. ||, is summarized in table | listing the 
longest series published until now [^,|6|,[7|,p|-p^1 f° r specific or generic values of N, on the sc 
and on the bcc lattices. 

This work is part of a sequence devoted to the extension and the analysis of HT series to 
order (3 21 for the iV-vector model on d- dimensional bipartite lattices. The case of the square 
lattice has already been discussed in Ref. [[J. We have chosen to use the (vertex renor- 
malized) linked cluster expansion (LCE) technique [|l3| , pOH 22]] because we have developed 
algorithms which make it equally efficient in a wide range of space dimensionalities indepen- 
dently of the nature of the site variables, whereas most other methods give best performances 
only in very specific situations, such as two-dimensional or low coordination number lattices 
and for sufficiently simple interactions. Having thus avoided the limitations of previous 
work, we have been able to produce extensive tables of series expansion coefficients given as 
explicit functions of the spin dimensionality N, which summarize in a convenient format a 
large body of information for an infinite set of universality classes. These tables, which we 
consider as the main result of our work, are reported in the appendices. 



We have built upon Ref. [13[] where some algorithms for the automatic LCE calculations 



have been introduced, and HT expansions of Xi °f 1^2 and of the second field derivative of 



the susceptibility Xi f° r the iV-vector model on the sc lattices have been tabulated |]T3| . [Td 
up to /3 1A for dimensions d = 2,3, 4. This calculation has been generalized to the class of 
d-dimensional bipartite lattices, in particular to the (hyper)sc and (hyper)bcc lattices, for 
which we have striven to design faster and more efficient algorithms and have introduced 
some innovations dramatically simplifying established computational schemes. In particular, 
by taking full advantage of the structural properties of these lattices, we have significantly 
reduced the fast growth of the combinatorial complexity with the order of expansion which 
had until now been the main obstacle to the extension of the series. In fact, it is mainly 
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our effort on the algorithmic side that has made progress possible, even in comparison with 
more recent work [23j using the same hardware resources. Moreover a considerable extension 



of our calculations is still feasible (and it is presently ongoing since we are far from 
our computational limits. Our calculation used an ordinary IBM RISC 6000/580H power 
station with 128 Mb memory capacity and 4 Gb of disk storage. Typical running times in 
the three-dimensional case were a few hours. In order to give a rough idea of the size of the 
computation it is sufficient to mention that over 2 x 10 7 graphs have to be generated and 
evaluated to complete the expansion of % an d ^2 through /3 21 . This should be compared 
with the corresponding figure: 1.1 x 10 4 , of the 0(/3 lA ) computation in Ref. |T3[ or with the 



figure 5 x 10 4 occurring in the recent analogous computation of Ref. |23| for the face centered 
cubic lattice to order /3 13 . Approximately 3 x 10 6 graphs contribute to the computation of 
X4 at order 17. 

We are confident that our results are correct, not only because our codes have passed 
numerous direct and indirect internal tests, but also because N and d enter in the whole 
computational procedure as parameters, so that a good general verification is achieved if 
our expansion coefficients, when specialized to iV = 0,1,2,3 and 00, agree with the (more or 
less) long series already available in various dimensions. Further details on the comparison 
with the available series, in particular in the limits N — > and iV — > 00, can be found in 
our paper devoted to the two-dimensional N— vector model 0. 

Note that, strictly speaking, the iV— vector model is defined only for positive integer 
spin dimensionality N. There are, however, infinitely many "analytic interpolations" in 
the variable N of the HT coefficients and, as a consequence, of all physical quantities. 
We have performed the "natural" analytic interpolation of the HT coefficients as rational 
functions of N, which coincides with that used in the 1/N expansion as well as in the usual 
Renormalization Group (RG) treatments and is unique in the sense of the Carlson theorem 



It is also worth emphasizing that the LCE technique can be readily adapted to produce 
HT expansions for the very general class of models (which include the O(N) symmetric 
P{0 2 ) lattice boson field theories), described by the partition function 

Z = J ' ILdn($)exp\p & ■ <Ps], (!) 

(r,s) 

where $ is a iV-component vector and dfi(0f) is the appropriate single spin measure. If 
we choose the form d[i{(pf) = 5{(pj — l)d<pi for the single spin measure, ([I]) reduces to the 
partition function of the iV-vector model. Also a broad class of other models of interest in 
Statistical Mechanics including the general spin S Ising model, the Blume-Capel model, the 
double Gaussian model, etc. can be represented in this form. The HT series for some of 



these models have been extended and will be discussed elsewhere p4|| . A wider discussion 
of ([1]) as a full theoretical laboratory for the study of scalar isovector lattice field theories 
as well as a detailed account of the graph theoretical and the algorithmic part of our work 
will also be presented elsewhere |24[ . 



The general interest in a direct determination of the critical properties of the classical 
lattice spin models with increasing reliability is clear. Other good general motivations for 
such a laborious calculation as a long HT series expansion include more accurate tests of 
the validity both of the assumption of universality, on which the Renormalization Group 
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(RG) approach to critical phenomena is based, and of the various approximation procedures 
required to estimate universal critical parameters by field theoretic methods. In fact, waiting 
for rigorous arguments to come, the only crucial tests [ES] of the validity of the Borel 
resummed e = 4 — d expansions [[26|-|28| or of the perturbative expansions at fixed dimension 
(FD) p6| , E8| -pl[1 , for which the O(N) model served as a paradigm, are presently limited to 
a careful comparison with experimental data or with other numerical data. Actually, this 
comparison mainly concerns numerical data of different origins, both because experiments 
are difficult and the experimental results, when available, are much less precise [with the very 
remarkable exception of the exponent v in the N = 2 case pifl ] and because experimental 
representatives are known only for a few universality classes [Q. Actually for N > 3 the 
physical, but not the numerical, interest of the model is somewhat lowered by the observation 
that the O(N) symmetric fixed point appears to be unstable within the e-expansion [33 



It should also be observed that, in spite of steady progress [jSSfl , the stochastic algorithms 
do not yet seem ready to completely supersede HT series in the study of models where the 
site variables have many components and/or the space dimensionality is large, or in the 
computation of multispin correlation functions. More generally HT series remain valuable 
subjects of independent study and sources of auxiliary information for other kinds of numer- 
ical calculations. Therefore we have also reported in the tables our estimates of nonuniversal 
critical parameters like the inverse critical temperatures. The computation of the critical 
amplitudes and of their universal combinations [36| will be discussed elsewhere [^4 . 



The paper is organized as follows: In section 2 we present our notation and define the 
quantities that we shall study. The analysis of the series is presented in section 3 along with 
a comparison to some previous analyses, to some results obtained by stochastic methods 
and to the RG results, both by the e-expansion and the FD perturbative techniques. In the 
appendices we have reported the closed form expressions for the HT series coefficients of \ 
and fi2 as functions of the spin dimensionality N and their evaluation for N = [the SAW 
model], N — 1 [the Ising spin 1/2 model], N = 2 [the XY model] and N = 3 [the classical 
Heisenberg model]. The present tabulation extends significantly and supersedes the one to 
order /3 1A in Ref. |14| which, unfortunately, contains a few misprints. 



In a forthcoming paper [EJJ] we will present an analysis of the series for the free energy 
and for its fourth field derivative \i- 



II. DEFINITIONS AND NOTATIONS 

For convenience of the reader we list here our definitions and notations. As the Hamil- 
tonian H of the A-vector model we shall take: 

H{v} = ~ £ v(x)-v(x'). (2) 

(x,x') 

where v (x) is a A-component classical spin of unit length at the lattice site x, and the 
sum extends to all nearest neighbor pairs of sites. 
The susceptibility is defined as 

oo 

X (N, (3) = 5>(0) ■ v(x)) c = 1 + E <*r(NW (3) 

x r=l 
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where (v(0) ■ v(x)) c is the connected correlation function between the spin at the origin 
and the spin at the site x. 

The second moment of the correlation function is defined as 

oo 

to{N,0) = £f>(0) ■ v(x)) c = J2 s rW r - (4) 

x r=l 

In terms of \ an d 1^2 we define the second moment correlation length £ by 



III. ANALYSIS OF THE SERIES 

Let us now turn to a discussion of our updated estimates for the critical temperatures and 
the critical exponents 7 and v in the N = 2, 3, 4 cases where our new series are significantly 
longer (up to 10 more terms) than those previously available. We shall also give some 
comments on the cases N = and 1 in which our extension is more modest and for N > 4 
where only a few numerical results are available. 

It has become clear from a long experience, mainly gained from the analysis of the Ising 
model HT expansions f6,37-41|, that in order to achieve a substantial improvement in the 



precision of the estimates of the critical parameters from the analysis of extended series 
one should properly allow for the expected nonanalytic corrections [§^] [usually also called 
confluent corrections] to the leading power law behavior of thermodynamic quantities near 
a critical point. For instance, we recall that, if we set r = 1 — /3//3 c , the susceptibility is 
expected to behave, in the vicinity of the critical point /3 C , as 

X {N, (3) ~ C x {N)t~^ (x + a x (N)r 6 ^ + a' x (N)r 26 ^ + ... + e x (N)r + e' x (N)r 2 + ...) (6) 

when r I 0. Not only the critical exponent 7(A), but also the leading confluent correc- 
tion exponent 0(N) is universal (for each N). On the other hand, the critical amplitudes 
C X (N), a x (N), a'(N), e x (N), etc. are expected to depend smoothly on the parameters of 
the Hamiltonian, i.e. they are nonuniversal. Similar considerations apply to £ (and to the 
other singular quantities) which, however, contains a different critical exponent and different 
critical amplitudes C^(N),a^(N), etc., but the same leading confluent exponent 6(N). It is 
also known that 6(N) ~ 0.5 for small values of N and 9(N) = 1 + 0(1/ N) for large N 

II- 

As experience has indicated, the established ratio extrapolation and Pade approximant 
(PA) methods are generally inadequate to the difficult numerical problem of determining 
simultaneously (3 C , the critical exponent and the leading confluent exponent in Eq. (^j), 
a task which essentially amounts to an intrinsically unstable double exponential fit. It is 
considered appropriate, then, to resort to the inhomogeneous differential approximants (DA) 



method |44j| , a generalization of the PA method, which, in principle, can be better suited to 
represent functions behaving like <f>i(x)(x — Xq)^ 1 + <p2 (x) near a singular point xq, where 
4>i(x) is a regular function of x and 02 (x) may contain a (confluent) singularity of strength 
smaller than 7. 
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A. Unbiased analysis 



To begin with, we have performed a series analysis by DA's, essentially following the pro- 
tocol suggested in Ref. [E_5] which is unbiased for confluent singularities. For each N, we have 
computed /3 C {N) and 7(A) by first and second order DA's built in terms of the susceptibility 
series and have then used this estimate of f3 c {N) to bias the determination of u(N) from the 
series for the square of the (second moment) correlation length £ 2 . Completely consistent 
results are also obtained, in general, by the method of critical point renormalization [03 



Also the specific heat exponent a(N) can be estimated by examining the behavior of x a t 
(3 = —/3 C (N), where a weak antiferromagnetic singularity is expected for bipartite lattices. 



As shown in Ref. ||46 |, having set f — 1 + /3//3 c , one has 

X (N, (5) ~ c{N) + ~a(N)f l ~ a{N) + b{N)f + ... (7) 



as f I 0. We shall however present this study in a forthcoming paper [[24j in order to 
jointly discuss also the results of the analysis of the free energy. 

The results of the present unbiased analysis do not significantly modify those obtained 
in the similar preliminary study || with series 0(j3 19 ). They are reported in table || for 
N < 3 and in table |TJ for N > 4 and compared with some of the most accurate recently 
published estimates P,p8|-p^.[r?|-|58| by various other methods, in particular by RG pertur- 



bative methods. Only for the physically most interesting cases, namely for N = 0,1,2,3, 



elaborate Borel resummed estimates are available for both the fifth order e-expansion [[2q -|28 
and either the six loop or the seven loop FD expansion |3Ti 



The scope of the seven loop FD computation of Ref. |3(J is however slightly limited 
by the present uncertainty in the value of the renormalized coupling constant g{N) used 



in the calculation. Therefore, in Ref. |ftO| , the exponent values have been conveniently 
expressed as the sum of the central estimate corresponding to this approximate value and 
of a small deviation proportional to the difference between g (N) and the true renormalized 
coupling g*(N). For simplicity, we have reported in table H t ne central estimate and have 
allowed for the contribution of the possible deviation only by doubling the expected error 
of the summation procedure: this roughly amounts to assume optimistically an uncertainty 
in the value of the renormalized coupling of a few parts per thousand. For N > 3, no 
estimates of the exponents by the e— expansion method have been published, while only 
very recently an extensive computation by the six loop FD expansion method has appeared 



| 31"| . Unfortunately, no estimates of error for the exponents are given in Ref. plf , but, in 



analogy with the small N case, we would reasonably expect relative errors of the order of 
1%. 

From our unbiased analysis of the sc lattice series we obtain exponent estimates essen- 
tially consistent, within their errors, with the available e— expansion results, but, in general, 
slightly larger (up to ~ 1%, or even more for intermediate values of N) than the FD expan- 
sion results. 

In the case of the bcc lattice, however, our unbiased estimates are also completely com- 
patible with the sixth order p6| , |29[| FD perturbative results or with the most recent seventh 



order |30| results. A similar situation has already been encountered in a previous very ac- 
curate unbiased analysis of the N = 1 case |j45| . We see this simply as a confirmation that 



the series for lattices with larger coordination number have a faster convergence rate [or, in 
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other words, a greater "effective length" H5]] and also as an indication that the simplest 



unbiased DA's might be only partially able to describe the confluent singularities. Therefore 
the larger discrepancies of the sc lattice results should not be interpreted as indicative of 
universality violations, but rather as a warning that the systematic errors of our analysis 
due to the finite length of our series and to the confluent singularities are quite likely to be 
underestimated if we evaluate the uncertainties solely from the scatter of the approximant 
values obtained using a sufficiently large number of series coefficients. We should also stress 
that, for N > 4, the sequence of DA estimates for the critical temperature or the exponents 
which use an increasing number of coefficients, show evident residual trends which indicate 
the presence of important confluent corrections to scaling, so that some "reasonable" extrap- 
olation of the results becomes necessary. Since this inevitably involves some assumption on 
the confluent exponent and therefore introduces some biasing, it will be more appropriately 
dealt with in the next paragraph. In order to distinguish clearly the effects of the various 
assumptions, we have chosen not to perform any further extrapolation in our "unbiased 
estimates" reported in tables || and [III], although it is clear that neglecting residual trends 
is a source of sizable systematic error. 

Even within these limitations, we have improved the precision of the values of the critical 
parameters as obtained from HT series by unbiased methods, and, so far, have not inferred 
from this analysis any indication of a serious inconsistency with the estimates from RG. 

In conclusion, within the unbiased approach to series analysis, the influence of the con- 
fluent singularities can be assessed more accurately and the results of the analysis can be 
better reconciled with the estimates from the RG methods probably only by computing still 
longer series, as it has been already recognized for the Ising model P, |38| , p9| . On the other 
hand, if we are ready to assume universality from the beginning, then, as suggested by work 
on the bcc lattice in the N = 1 case whose results are reported in table [TT], a more 

accurate determination of the critical exponents could be obtained, even without further 
extending the series, by a study of appropriately built families of models depending on some 
continuous parameter, for each universality class. The idea is, essentially, to minimize or 
suppress the amplitude of the dominant confluent correction to scaling in Eq. (|6]) by tak- 
ing advantage of its nonuniversality, namely of its continuous dependence on the parameter 
entering in the Hamiltonian. Using our LCE computation, it is now possible to implement 
easily this procedure for any N and on two different lattices, thereby corroborating its re- 
liability. These developments of our study and a detailed analysis of other features of the 
series including estimates of the universal ratios of the leading confluent amplitudes for x 
and £ 2 will be presented elsewhere p4 . 



B. Biased analysis 

We have also analyzed our series by various biased methods, in particular by using 
properly designed first order inhomogeneous DA's in which both (3 C {N) and the correction 
to scaling exponent 0(N) have been fixed , or by second order inhomogeneous DA's in which 
only 9(N) 0J5| has been fixed. 

In order to provide the additional information needed in these approaches, for N < 3, 
we have assumed that the exponents 0(N) take the values predicted by the FD perturbative 
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9(0) = 0.470(25), 9(1) = 0.498(20), 9(2) = 0.522(18), 9(3) = 0.550(16). (8) 



For iV > 3, the first FD perturbative estimates at six loop order have only been ob- 
tained very recently by A. I. Sokolov and have been kindly communicated to us before 
publication 

0(4) = 0.578(10), 9(6) = 0.626(10), 9(8) = 0.670(10), 9(10) = 0.707(10), 9(12) = 0.737(10). (9) 

We have therefore been able to revise the biased analysis presented in the first preprint 
version of this paper, where, for lack of a better choice, we had used values for 9(N) (N > 3) 
(larger than those in Eq. (|9])), obtained by a reasonable, but unwarranted interpolation 
method. In some cases, in particular for N = in Refs. P, f4"T| , |6"Tf and for iV = 1 in Refs. 
|7|j62]|, numerical work has suggested slightly different confluent exponents, which might be 



more accurate or, to some extent, also provide an "effective" description of higher confluent 
corrections. We shall return later to this point, but we should mention that, at the level 
of precision of the following calculations, the precise values of the exponents 9(N) do not 
matter too much, since we have observed that our biased estimates of the leading critical 
exponents remain practically stable (within their errors) under variations of the confluent 
exponents up to ~ 5 — 10%. Of course other quantities, such as the confluent amplitudes, 
are very sensitive to the values of 9(N). 

Let us now sketch the first biased method. If an accurate estimate for f3 c is available, we 
can formulate a quite simple procedure, biased with both /3 C and 9. From the asymptotic 
formula (§) we get 

7 (iV, P) ee (P C (N) - ^ dl0 ^(N,(3) = ^ {N) + p(Nm{N) _ ^(N) + om(N) _ p )9 W) 

(10) 

where p(N) = -9(N)a x (N). 

We can approximate the quantity j(N, (3) by the particular class of first order inhomo- 
geneous DA's defined as the solutions of the equation 

QM [((3 C (N) + 9(N)y(f3)} + R n ((3) = (11) 

Here Q m (P) and R n (/3) are polynomials of degrees m and n respectively, calculated, as 
usual, from the known series expansion of ^(N,^). As a result the exponent 7 is simply 
estimated as 

(m -Rn((3 c (N)) 



9(N)Q m (Pc(N)) 

and the amplitude p(N) of the subleading term in Eq. ( |TDD is given by the formula 

n(m y(0)-7(AQm,n f^ N) D(t)m,ndt 

P[ )m ' n {3 c (Nf N) J° ((3 c (N)-t) 1+9 ^ [ ] 

where 
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lUtl _ RMN)) 

[)m ' n Qm(t) Qm(f3c(N)) 

We consider only almost diagonal approximants, namely those with \m — n\ < 4. This 
procedure, which might be seen as the simplest, although certainly not the most general, 
DA extension of the biased PA method introduced in Refs. p3| , |64| , works accurately on 
model series having the analytic structure expected for j(N, f$) in the vicinity of f3 c (N). 
A similar procedure can be applied to £ 2 (N,/3)/fl in order to compute the exponent u(N). 
In order to give an idea of the results, let us for example, consider the A = 1 sc lattice 
series. Assuming (3 C = 0.2216544(3) as suggested in Ref. |49j and 9 = 0.498(20), we estimate 



7 = 1.2388(10) and v = 0.6315(8). With the same value of 9, in the bcc lattice case, assuming 
(3 C = 0.157373(2) as suggested in Ref. 0, we obtain 7 = 1.2384(6) and v = 0.6308(5). In 
general, when the value of (3 C is not known accurately, it would be more appropriate to 
present the results in the form of a linear relationship between the critical exponent and /3 C , 
at a given fixed value of 9. We can, however, use this method also to determine the value 
of (3 C by fixing only 9 and looking for the (generally small) range of values of (3 C for which 
the uncertainty of the exponent 7 is minimal. The estimates obtained in this way are in 
general completely consistent with those from the analogous biased PA procedure of Refs. 
|63| , |64]| . The values of (3 C so obtained are also, in most cases, consistent with those obtained 



from the unbiased improved ratio methods of Ref. |38[ after extrapolating the sequences of 
results linearly in l/r 1+9 ^ (where r is the number of series coefficients used), and with the 
estimates obtained by similarly extrapolating the results from unbiased DA's. Analogous 
considerations apply to our exponent estimates, which, however, have to be compared to 
the results obtained from improved ratio methods by extrapolating linearly in l/r e ^ N \ 

A second method, modelled after Refs. [[7|,pU|, which is biased only with the value of the 
confluent exponent, may be described as follows: for each value of N, we have considered 
the approximants derived from inhomogeneous second order differential equations with the 
structure 

(&(A0 - (3) 2 ((3 C (N) + P)QiW^jji + Wo(N) - P)Pm(P)^ + RniPW) + T s ((3) = 

where Qi(fi), P m (j3), R n {(3) and T s {(3) are polynomials of degrees /, m, n and s in the variable 
(3. By this choice the DA's are biased to be singular at (3 = /3 C (N) and at (3 — —f3 c (N). 
We restrict ourselves to almost diagonal DA's (namely those with / + 3~m + l~n and 
s < 4), which use at least 19 series coefficients. For each DA, we adjust (3 C (N) in a small 
range around the values indicated by the unbiased analysis of the previous section until the 
correction to scaling exponent 9(N) reaches precisely the central value indicated in Eqs. 
(!);©■ The corresponding inverse critical temperature and exponents are then taken as 
the best estimates of these quantities. It should be noticed that, within this approach, the 
values of (3 C for x an d £ 2 cannot be forced to be equal, but the differences are generally not 
much larger than the errors. In this approach (like in the previous one) the evaluation of 
the errors is based, as usual, on the spread among the approximant values and includes a 
generous allowance for the uncertainty in the biased value of 9(N). Let us finally mention 
that we also have preliminarily tested the reliability of this procedure on various model series 
built in such a way to reproduce the main expected features of x an d £ 2 - 
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An important remark concerns the values of the amplitudes of the confluent corrections 
in x an d £ 2 , which can be obtained from Eq. (fL2|). We have observed that, in the bcc 
lattice case, they are negative for N < 2, positive and increasing for N > 2. A quite 
similar behavior of these amplitudes is observed in the sc lattice case, where the change 
in sign occurs for N ~ 3. For 2 < iV < 4, the amplitudes are relatively small for both 
lattices, but they become important for N > 4. This remark, which will be illustrated in 
full detail elsewhere pi |, is also consistent with the fact that our biased estimates of the 
critical exponents are increasing functions of the confluent exponents for N < 2, while they 
are decreasing functions for N > 4. For 2 < N < 4 some approximants give increasing 
functions and others give decreasing functions. Completely consistent indications on the 
behavior of the confluent amplitudes can also be inferred from the features of the improved 
ratio plots J38], or by studying, as function of N, the difference among the unbiased DA 



estimates of (3 C obtained from x, /i 2 or £ 2 , as suggested in Ref. [f39fl . In the N = 1 case, our 



remark agrees with the results of earlier studies |7| j3^| , |Ho1^T| on the sc, bcc and fee lattices, 
which established that the sign of the leading confluent amplitudes in \ an d £ 2 is negative. 

Some features of the biased procedures we have adopted, may appear questionable or 
may call for further improvement. We have always tried to minimize the possible defects by 
forming our final estimates with the, generally quite compatible, indications coming from all 
the various available biased methods and not only from the two biased DA procedures we 
have described above. Moreover we have indicated very conservative error bars. It is however 
reasonable to expect that the accuracy of the biased estimates depends not only on whether 
the present series are long enough to provide the essential information on the subdominant 
singularities, but also on whether most of the corrections to scaling can actually be well 
described by the first nonanalytic term in the asymptotic formula (^|). We believe, however, 
that the final results are at least consistent and suggestive. First of all, our analysis confirms 
that the total (statistical + systematical) errors of the previous unbiased approach are larger 
than we have indicated. For N < 1 the biased exponent estimates differ only slightly from 
the unbiased ones and in such a way to improve the agreement with the most accurate 
FD perturbative estimates, with stochastic simulations and with experimental results. For 
N = 2, 3 the agreement with the RG, in the FD perturbative approach, is, perhaps, less 
convincing. For N > 4 our biased exponent estimates are systematically larger (up to ~ 4%) 
than the FD six loop perturbative values [BT]]. This discrepancy is parallel to and, in our 



opinion, is strictly related to the fast increase of the confluent amplitudes of x an d £ 2 in the 
same range of N. Therefore the problem seems to be of a purely numerical nature. Of course, 
we cannot prove that, in these circumstances, our biased analyses can allow for the confluent 
singularities better than the present Pade-Borel resummation techniques do within the FD 
perturbative approach. We can perhaps only argue that, at least for large N, our results 
might be more accurate by comparing them with the results of the 1/N expansion of the 
exponents. Indeed for N > 14, the 1/N estimates, which seem to be rather well converged 
and therefore reasonably accurate, seem to approach our biased bcc estimates faster than the 
FD values. Let us finally add that these numerical problems with the perturbative estimates 
of the exponents might also be partly due to possible residual imprecisions in the values of 



the renormalized coupling constants entering in the FD perturbative calculation [Q. We 
shall return to this point when new estimates of the coupling constant will be available from 
HT series lEl. 
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C. Final comments 



Let us now finally add some comments on the results of the series analysis which are 
presented in the tables. 

In the N = case, the SAW model, we have not attempted either to report in table 
1 1] or to cite in the references even only a representative sample of the large amount of 
numerical work accumulated over the years, which fortunately has been reviewed recently in 



the very extensive and valuable new treatise |66j and in Ref. |47J devoted to a high precision 
stochastic study. 

On the sc lattice, our series for x and /12 are not yet longer than those of Ref. || and of 
Ref. respectively, but we have taken advantage of the two additional published expansion 
coefficients of \ in order to test the stability of our biased estimates. We recall that the 
previous HT analysis by unbiased DA's of the sc lattice series to order (3 21 performed in 
Ref. § produced the estimates /3f = 0.213496(4), 7 = 1.161(1) and v = 0.592(4), which 
are all completely consistent with the results of the later analysis of Ref. |4]] and with our 
own unbiased analysis, but slightly larger than our biased estimates. Very recently the 
significantly lower estimates f3 s c c = 0.213492(1) and 7 = 1.1575(6) have been obtained by a 
stochastic method sampling SAWs which extend up to 4 x 10 4 steps |48j . Using this estimate 
of j3 s c c in the first DA method introduced in the previous section, we get 7 = 1.1589(8); in 
order to obtain the lower value 7 ~ 1.1575 a value (3 S C C ~ 0.213488 would be required, which 
is somewhat below the range suggested by the presently available HT series on the sc lattice, 
even allowing for the confluent corrections. 

In the bcc lattice case we have computed five new coefficients for \ and fi2 beyond those 
reported in Ref. [[J. This makes it worth computing new estimates for the exponents. We 
recall that the analysis in Ref. of the 0(f3 m ) bcc lattice series available until now yielded 
the values (3 b c cc = 0.153137(10), 7 = 1.162(2) and v = 0.592(2), which are less precise, but 
compatible with our new unbiased estimates and somewhat larger than our corresponding 
biased estimates. On the other hand, our biased estimate of 7 for the bcc lattice agrees more 



closely with the new stochastic estimate ||48|| . We should also stress that, for both the sc and 
the bcc lattices, our biased estimates of v have now come very close to the RG estimates 
and to the experimental value v = 0.586(4) reported in Ref. |)7 . 



The value of #(0) is still controversial | ]61[| . For example, the study j|] of the long HT 
series available on the sc lattices suggests #(0) ~ 1., while an extensive MC study |47 



on the same lattice rather indicates an effective exponent #(0) ~ 0.56(3). Assuming this 
last value in our computation instead of the one in Eq. ([8]), would only shift the biased 
estimates of 7 and u, in the bcc lattice case, from 1.1582(8) to 1.1587(8) and from 0.5879(6) 
to 0.5883(6) respectively. Similarly in the sc lattice case the estimate of 7 (using all 23 
available coefficients) would change from 1.1594(8) to 1.1597(8) and that of v from 0.5878(6) 
to 0.5894(8). 

A final remark is that from our extended series for x on the bcc lattice we can derive 
|5T| the new rigorous and stricter inequality /^ cc (0) > exp(— ^-/na2i(0)) ~ 0.148582.., which 
slightly improves the previous bound obtained from ai 6 (0) and quoted in Ref. [Q. 

In the N = 1 case, the Ising spin 1/2 model, the relevant numerical studies are even more 
numerous than for the SAW model, so that we can only address the reader to the recent 
extensive review which in Ref. @] complements stochastic computations of unprecedented 
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accuracy on the sc lattice. On this lattice, we have extended by two terms the \ series and 
by six terms the /i 2 series [ 1 1 so that it is worthwhile at least to update the estimate of v. 
In order to compute the exponents reported in table [Tl] from the sc series, we have assumed 
9(1) = 0.498(20) and simply taken the extremely accurate value /3f (1) = 0.2216544(3) 
obtained in Ref. |0J and therefore indicated within parentheses. Our bcc lattice series for 
X and /i 2 are not yet longer than the series of Refs. In this case, we have assumed 

the same value of 0(1) and taken j3^ cc (l) = 0.157373(2) from Ref. M (similarly indicating it 
within parentheses). Also for the Ising model, values of the confluent exponent slightly larger 
than the one we have assumed, such as #(1) = 0.52(3) in Ref. J7J and #(1) = 0.54(3) in Refs. 
39| , |62[| have been reported. Taking the largest of these values in our biased computation 



would only change, in the bcc lattice case, the central estimate of 7 from 1.2384 to 1.2387 
and would shift that of v from 0.6308 to 0.6311. In the sc lattice case the central value of 7 
would be shifted from 1.2388 to 1.2392 and that of v from 0.6315 to 0.6321. 

In the N = 2 case, the XY model, we have computed four more terms of the x an d At 2 
series in the sc lattice case. Notice that the last two coefficients of the previously published 
0(f3 17 ) series contained tiny numerical errors (inconsequential for the analysis in Ref. |fT2| ) 
which are corrected by our new computation. In the case of the bcc lattice, our extension 
of the series for x and yU 2 amounts to nine terms and gives a greater significance to the new 
exponent estimates. 

In the N = 3 case, the classical Heisenberg model, we have extended by seven terms the 
series for x an d ^2 on the sc lattice. In the case of the bcc lattice we have extended the 
series by ten terms. 

For all iV > 3, only series up to 0(/3 9 ) were available until now on the bcc lattice and 
therefore our extension amounts to twelve terms. On the sc lattice we have computed seven 
additional series coefficients for x an d /-t 2 . Let us finally note that, on the sc lattice for 
N = 2, 3 and 4, the estimates of j3 c indicated by the simulation of Ref. ]54j] are only slightly 
smaller than ours. Using these values in our biased DA method, we would obtain 7 ~ 1.322, 
7 ~ 1.402 and 7 ~ 1.476, for N = 2, 3 and 4, respectively. For the exponent v the agreement 
would be somewhat closer. 



IV. CONCLUSIONS 

In conclusion, we have produced new HT expansions through order /3 21 of the suscepti- 
bility and of the second correlation moment for the classical iV-vector model with general 
N, on the sc and the bcc lattices. This rich material has been conveniently tabulated in the 
appendices in order to offer an easy opportunity for further study. 

As a first application of our results, we have updated the direct estimates of the critical 
parameters of the N- vector model with a considerable improvement in accuracy over previous 
analyses and, for all values of N, we have confirmed a generally good agreement with the 
most precise calculations by current approximate RG methods. 
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TABLES 



TABLE I. Longest published HT expansions for the N- vector model on the simple cubic and 
the bcc lattice before this work, 

Quantities expanded Parameters Maximal order Ref. 

sc lattice 

X N = 23 | 

X ,H2 N = 21 § 

X N = l 19 

N = l 15 



M2 



X,V2 N = 2 17 [12 



X,V2 anyiV 14 MM\ 



bcc lattice 

X,M2 iV = 16 H 

X^2 N = l 21 |,0 

X^2 N = 2 12 § 

X N = 3 11 @ 



X,/i 2 anyiV 9 @ 
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TABLE II. A summary of the estimates of the critical parameters for < N < 3 


N 


Method and Ref. 


Pc 


7 


V 





HTE sc unbiased [§ 


0.2134987(10) 


1.16193(10) 






HTE sc unbiased 


0.213497(6) 


1.161(2) 


0.592(2) 




HTE sc 0-biased 


0.213493(3) 


1.1594(8) 


0.5878(6) 




MonteCarlo sc Mjj 


0.2134969(10) 




0.5877(6) 




MonteCarlo sc J4q] 


0.213492(1) 


1.1575(6) 






HTE bcc unbiased 


0.153131(2) 


1.1612(8) 


0.591(2) 




HTE bcc 0-biased 


0.153128(3) 


1.1582(8) 


0.5879(6) 




R.G. FD perturb. Q 




1.1569(8) 


0.5872(8) 




R.G. e— expansion |2o| 




1.157(3) 


0.5880(15) 


1 


HTE sc unbiased 


0.221663(9) 


1.244(3) 


0.634(2) 



HTE sc 0-biased 
MonteCarlo sc [ j£ ] 
MonteCarlo sc j E \ 
HTE bcc 
HTE bcc ||r 
HTE bcc unbiased 
HTE bcc 0-biased 
R.G. FD perturb, j 
R.G. e— expansion 1 



(0.2216544(3)) 
0.2216595(26) 
0.2216544(3) 



0.157379(2) 
(0.157373(2)) 



1.2388(10) 

1.237(2) 

1.2395(4) 

1.237(2) 

1.243(2) 

1.2384(6) 

1.2378(12) 

1.2390(25) 



0.6315(8) 

0.6289(8) 

0.6301(8) 

0.632(1) 

0.6300(15) 

0.634(2) 

0.6308(5) 

0.6301(10) 

0.6310(15) 



Experiment [ p2[ 
HTE sc unbiased 
HTE sc 0-biased 
MonteCarlo sc [ y ] 
MonteCarlo sc [ % \ 
MonteCarlo sc j oA j 
HTE bcc unbiased 
HTE bcc 9-biascd 
HTE fee (X5| 
R.G. FD perturb. [ 
R.G. e— expansion 



0.45419(3) 

0.45419(3) 

0.45420(2) 

0.4542(1) 

0.454165(4) 

0.320428(3) 

0.320427(3) 

0.2075(1) 



1.327(4) 

1.325(3) 

1.308(16) 

1.316(5) 

1.319(2) 

1.322(3) 

1.322(3) 

1.323(15) 

1.318(2) 

1.315(7) 



0.6705(6) 

0.677(3) 

0.675(2) 

0.662(7) 

0.670(7) 

0.672(1) 

0.674(2) 

0.674(2) 

0.670(7) 

0.6715(15) 

0.671(5) 



HTE sc unbiased 
HTE sc 0-biased 
MonteCarlo sc fea ] 
MonteCarlo sc p4] 
HTE bcc unbiased 
HTE bcc 6»-biased 
MonteCarlo bcc |§ 
HTE fee jl6) 
R.G. FD perturb. [M 
R.G. e— expansion p6| 



0.69303(3) 

0.69305(4) 

0.693035(37) 

0.693002(12) 

0.486805(4) 

0.486820(4) 

0.486798(12) 

0.3149(6) 



1.404(4) 

1.406(3) 

1.3896(70) 

1.399(2) 

1.396(3) 

1.402(3) 

1.385(10) 

1.40(3) 

1.3926(26) 

1.39(1) 



0.715(3) 

0.716(2) 

0.7036(23) 

0.7128(14) 

0.711(2) 

0.714(2) 

0.7059(37) 

0.72(1) 

0.7096(16) 

0.710(7) 



18 



TABLE III. A summary of the estimates of the critical parameters for 4 < A < 12 



N 


Method and Ref. 




/3 C 


7 


V 


4 


HTE sc unbiased 




0.93589(6) 


1.474(4) 


0.750(3) 




HTE sc 0-biase< 






0.93600(4) 


1.491(4) 


0.759(3) 




MonteCarlo sc 






0.9360(1) 


1.477(18) 


0.7479(90) 




MonteCarlo sc 






0.935861(8) 


1.478(2) 


0.7525(10) 




HTE bcc unbiased 




0.65531(6) 


1.461(4) 


0.744(3) 




HTE bcc 0-biased 




0.65542(3) 


1.484(4) 


0.756(3) 




R.G. FD perturb. 


|3 




1.45(3) 


0.74(1) 




R.G. FD perturb. 


ifj 




1.449 


0.738 


6 


HTE sc unbiased 




1.42859(6) 


1.582(5) 


0.804(3) 




HTE sc 0-biased 




1.42895(6) 


1.614(5) 


0.821(3) 




HTE bcc unbiased 




0.99613(6) 


1.566(4) 


0.796(3) 




HTE bcc 9-biased 


1 — 1 


0.99644(4) 


1.608(4) 


0.819(3) 




R.G. FD perturb. 




1.556 


0.790 


8 


HTE sc unbiased 




1.9263(2) 


1.656(5) 


0.840(3) 




HTE sc 0-biased 




1.92705(7) 


1.701(4) 


0.864(3) 




HTE bcc unbiased 




1.33984(7) 


1.644(5) 


0.833(3) 




HTE bcc S-biased 


1 1 


1.34040(6) 


1.696(4) 


0.862(3) 




R.G. FD perturb. 






1.637 


0.830 




1/N expansion 


M 






1.6449 


0.8355 


10 


HTE sc unbiased 




2.4267(2) 


1.712(6) 


0.867(4) 




HTE sc 0-biased 




2.42792(8) 


1.763(4) 


0.894(4) 




HTE bcc unbiased 




1.68509(8) 


1.699(5) 


0.860(4) 




HTE bcc 0-biased 




r.uooou^f j 


1.761(4) 






R.G. FD perturb. 






1.697 


0.859 




1/N expansion 








1.7241 


0.8731 
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HTE sc unbiased 




2.9291(3) 


1.759(6) 


0.889(4) 




HTE sc 0-biased 




2.9304(1) 


1.812(5) 


0.916(4) 




HTE bcc unbiased 




2.03130(8) 


1.741(6) 


0.881(4) 




HTE bcc 0-biased 




2.03230(8) 


1.808(5) 


0.914(3) 




R.G. FD perturb. 






1.743 


0.881 




1/N expansion 


m 






1.7746 


0.8969 
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